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Abstract 



We study time dependence of various measures of entanglement (covariance entanglement coefficient, 
purity entanglement coefficient, normalized distance coefficient, entropic coefficients) between resonantly 
fS) , coupled modes of the electromagnetic field in ideal cavities with oscillating boundaries. Two types of 
cavities are considered: a three-dimensional cavity possessing eigenfrequencies W3 — 3uJi, whose wall 
oscillates at the frequency luw = 2uJi, and a one-dimensional (Fabry-Perot) cavity with an equidistant 
^ . spectrum ujn = nuji, when the distance between perfect mirrors oscillates at the frequencies loi and 
2ijJi. The behaviour of entanglement measures in these cases turns out to be completely different, 



. although all three coefficients demonstrate qualitatively similar time dependences in each case (except 
Q,^ ■ for some specific situations, where the covariance entanglement coefficient, based on traces of covariance 
' submatrices, seems to be essentially more sensitive to entanglement than other measures, which are based 
on determinants of covariance submatrices). Different initial states of the field are considered: vacuum, 
• squeezed vacuum, thermal, Fock, and even/odd coherent states. 



^ PACS: 42.50.Lc; 42.50.Dv; 03.65.-w 

Key words: Dynamical Casimir effect; Vibrating boundary; Parametric resonance; Coupled modes; En- 
Jrtanglement; Quantum purity; Entropy; Distance; Covariances; Fock states; Gaussian states; Even/odd 
Coherent states; Squeezed states; Thermal states 

cr. 

Introduction 

K . 

J-JDUring the past decade it was recognized that the concept of entanglement^ introduced by Schrodinger in 
. 5^935 [0, ^, is not only one of the most profound in quantum mechanics (as was shown in the same year by 
Einstein, Podolsky and Rosen in their famous paper albeit without using explicitly this word), but it is 
crucial for many promising new applications, such as quantum cryptography, quantum communication and 
teleportation, quantum computing, etc. This explain a burst of interest to various problems connected to 
this concept observed for the past few years. One of such problems is a search for quantitative measures of 
entanglement. 

In the most cases, the measures based on different kinds of entropies have been considered ||, H, ^, ||, ^ . 
For example, if the total system, consisting of parts 1 and 2, is decribed by means of the statistical operator 
po, then the entanglement measure is frequently expressed in terms of the total and "partial" entropies as 
the "index of correlation" 

Ic = Si + S2 ~ So, Sk = -Trfe pk Inpk, (1) 

where the reduced statistical operator is defined as, e.g., pi — Tr2po- For the total pure states formula (|^) 
is reduced to Ic ~ 2Si — 2S2- 

*on leave from Lebedev Physical Institute and Moscow Institute of Physics and Technology, Russia 
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However, despite of many advantages, the measures such as (|l|) are not very convenient from the practical 
point of view: in order to calculate them, one has to diagonalize the reduced statistical operators, and this is 
rather difficult problem in the generic case, especially for infinite-dimensional Hilbert spaces (corresponding 
to the so called "continuous variable systems"), except for a few simple special cases. Therefore, many 
authors looked for other measures, which could be calculated more easily. 

In our paper, we consider several families of the simplified measures. The first one is based on the notions 
of quantum purity /i = Trp^ or "linear entropy" = 1 — /i. Different measures containing these quantities 
were proposed in ||, ||, |lO[ . Measures based on the Hilbert-Schmidt distance between the given state 
and its "disentangled" counterpart were proposed in |3j. On the other hand, even simpler (although 
non-universal) measures of entanglement of continuous variable quantum systems, expressed in terms of the 
cross-covariances of the quadrature components or the annihilation/creation operators, have been introduced 
recently in jl4| These measures are discussed in Section 2. 

One of numerous possible applications of the entanglement measures is a compact quantitative charac- 
terization of the evolution of coupled quantum mechanical systems. We began these studies in Q, where 
two harmonic oscillators with constant frequencies but with the most general time-dependent resonance 
couplings were considered. The aim of the present paper is to compare different measures in the case when 
the field modes in a cavity are entangled due to the motion of its boundary (the physical reason of entangle- 
ment in this case is the Doppler effect). This case is reduced to the models of two or many oscillators with 
time- dependent frequencies and a specific time-dependent coupling (of the "coordinate-momentum" type). 

We consider two types of cavities, beginning with a three-dimensional cavity with accidental degeneracy 
of the spectrum (which happens, e.g., in cubical cavities), when only two modes can occur in resonance with 
an oscillating wall (Section 3). A one-dimensional (Fabry-Perot) cavity is considered in Section 4. In this case 
all modes are coupled due to the equidistance of the (unperturbed) spectrum of the field eigenfrequencies. 
It was discovered as far back as in [p^ that the field evolution in three- and one-dimensional cavities is 
qualitatively different. For example, in the 3D case the number of photons in the resonance modes grows 
with time exponentially, whereas in the ID case this growth is only linear. It was pointed out in that 
the growth of the number of photons in the ID cavity is slowed down due to a strong intcrmode interaction, 
which is equivalent in this case to entanglement. Now we are able to give a quantitative characterization of 
such an entanglement. The results of our study are discussed in Section 5. 

2 Purity, distance and covariance measures of entanglement 
2.1 Purity entanglement measure 

By analogy with definition (^, the "linear entropy of entanglement" can be defined as 



Such a definition seems reasonable if the total system is in a pure quantum state. Then Trp^ — 1 and 
Tri pf = Ti2 P2, so that C — 25f = 25f . As a matter of fact, only this case was considered in the earlier 
studies where measures of entanglement were identified with the linear entropy of the state of a 

subsystem or with some equivalent quantities, such as the purity itself, the "participation ratio" l/Trfc/5^, 
or the "Renyi entropy" S'-'^ = — In (Tr^ p^) . 

However, if the state of the total system is mixed, then definition leads to some unexpected con- 
sequences. Consider, for example, a generic Gaussian two-mode state described by means of the Wigner 
function (we assume h= 1 throughout the paper) 



£ = 5f + - = 1 + Tr f? - Tn pi - Tr2 pi 



(2) 



W(q) = |det(Q)ri/^xp [-1 (q- (q» Q-^ (q- (q)) 




(3) 



where q = [xi,pi, X2,P2), and the symmetrical 4x4 covariance matrix Q consists of 2 x 2 blocks 



S= lka/3| 



Qii Qi2 
S21 S22 



Qii — Qii, Q22 — Q22 



Q12 = Q21 



(4) 
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(a tilde over matrices means matrix transposition). The symmetrical real covariances are defined as 

qap = 2 (9^^ + qMa) = <iaqi3, ab = (aS) - (a) (6) (5) 

(in other words, a straight line over the product of two observables means the ordered centralized average 
value, whereas a wide tilde means the symmetrized centralized average value). Then 

fi^TTp^=J [W{ci)fdq/{27Tf = [det(2Q)]-i/2 . (6) 

For factorized (disentangled) states, Q12 = 0, therefore det Q — det Qn det Q22 and /i = /ii/i2, which 
results in the relations 

Cfact = 1 + - Ml - = (1 - Ml) (1 - M2) • 

Consequently, for mixed states (/i < 1), one can meet the situation when £ > in the absence of any 
entanglement, if /ii 7^ 1 and /i2 7^ 1- 

It seems better to use the difference C^, = C — Cfact = M ~ MiM2- But it tends to zero when — > 0. For 
this reason, we introduce the normalized purity entanglement coefficient 

L^l-f^. (7) 
M 

For the Gaussian states (p[) it can be expressed as 



detQ 



^ = ^ - V detg^dltQ.. - ' - V^^t - 2-2.-.^2-2n^) ■ (8) 

The second equality (where E stands for the unit matrix) is obtained with the aid of the known formula for 
the determinant of a block matrix |]l7|| 

det Q = det (Qn - Qi2Q22'22i) det Q22. 

In particular, for pure composite states (/i = 1) we have 

£ = l-/i? = l-M2 = J^(4-/:). (9) 

A measure of entanglement between two coupled modes resembling (^ was introduced in [ p^ (where it 
was named "group correlation coefficient"): 

^ =C(2-C). (10) 
detQiidetQ22 V ^ ^ ' 

In principle, the measures (||) and ( |lO| ) can be used for arbitrary (not only Gaussian) states, although 
sometimes they can give zero value even for truly entangled states (if the matrix of the second-order variances 
is factorized, but intermode correlations exist for higher-order moments). Also, instead of (j7|) one could use 
the following extension of formula (^0|) to arbitrary states: 

=2 -, /MiM2^^ 



IC' = l^i^^j ^^^^ 
2.2 Distance entanglement measure 

Another possibility to characterize entanglement is to use the Hilbert-Schmidt distance between the given 
state and different "disentangled" states. It was considered, e.g., in [Q, ^ ^ (analogous approach was 
developed in [|l9) to quantify the "degree of nonclassicality" of quantum states). In the entanglement 
measure was defined as Tr {p — pi® p-if ■ However, we prefer to normalize it by Trp^, in order that the 
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entanglement measure would not go to zero for highly mixed states. Thus we shall consider the following 
quantity: 



Tr [p-pi® P2) 
Trp^ 



l + ^--Tr(/5-[pi®p2]). 



(12) 



For any states p and R one has (the normalization factor corresponds here to the two-mode case) 

Tr(pi?) = j Wp(q)W^fl,(q)dq/(2^)2 . (13) 
For the Gaussian states (^ the integrals can be calculated with the aid of the known formula 

J exp (-qAq + bq) dq = [det(A/7r)]~^/^ exp QbA^^b^ . (14) 
The inverse matrix can be represented in the block form with the aid of the Frobenius formula 

(15) 

Q* = Q22- Q2iQn 2i2- 

Taking into account equations (||), (||), ^4), and (|5|), one can verify that the Wigner function of the 
factorizcd state pi (8" P2 is given by formula ( 3|) with the block-diagonal matrix 



Qii 


Q12 


-1 




Q21 


Q22 




-Q-^Q2iQii Q-^ 



Qii 
Q22 



where matrices Qn and Q22 are the same as in (||) (this is obvious from the physical point of view). Thus 
we arrive at the following expression for the 2-measure (it is equivalent, except for the normalizing factor 
/i^^, to that given in |13 , but it is written in more simple explicit form): 



Z = 1 



detQ 



det Qn det Q2 



/det(2Q) 
detQ, ' 



Q+Qd 



2Qii Q12 
Q21 2Q22 



(16) 



2.3 Covariance entanglement measures 

Other measures of entanglement have been introduced recently in They are expressed directly in 

terms of the cross-covariances of the quadrature components or the equivalent annihilation/creation operators 
as follows: 



y 



Tr(Qi2Q2i) 
TrQiiTrQ22 



1/2 



|ai4p + \aia2\ 



1/2 



{xiX2f + {piP2f + {xiP2f' + {PlX2f 
4S1S2 



1/2 



y 



2 {a\ai + 1/2) {a\a2 + 1/2^ 
2v/Tr(Qi2Q2i) _ ^2(|ai4p + |aTai|2) ^(3^)^ + (|^) V (a^) V (|^) 



where (we use properly normalized dimensionless quadrature variables) 

T- 1 1 



afc = (ife + ipk) 



£k = + § = § {^kXk + PkPk) 



fc = 1,2. 



(17) 

(18) 

(19) 
(20) 
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Since the coefficients (|8|) and are expressed in terms of traces of products of the off-diagonal blocks of 
the total covariance matrix Q, they are obviously invariant with respect to the rotations in the phase plane 
of each subsystem. (Another invariant quantity, namely the determinant of the off-diagonal blocks, det Q12, 
plays an important role for the problem of separability of continuous variable systems po| ). It can be shown 
that < < 3^ < 1. 

We would like to emphasize that the coefficients y and y are defined for any (not only Gaussian) 
quantum state. They are significantly simpler than other entanglement measures from the point of view of 
calculations (to calculate traces of matrices is much more easy than to calculate determinants, not speaking 
on calculating eigenvalues of density operators or matrices, which are necessary to obtain the entropic 
measures). A disadvantage of the coefficients y and y is that in some cases they are equal to zero even when 
the state is entangled, but the second-order moments of quadrature components are equal to zero. However, 
this does not happen for Gaussian and many other important quantum states. 



2.4 Entropic measures for Gaussian states 

In order to demonstrate how simple are expressions given in the preceding subsections, compared with the 
"standard" entropic measure (|l|), we give here the formula for the entropy of a generic Gaussian state. It is 
also determined by the covariance matrix, but in a more complicated way than the coefficients considered 
above. For an arbitrary iV-mode Gaussian state the entropy was found in different but equivalent forms in 



|2l| , |22|] and recently in [g3| . The most simple expression is 1 24 

^{kj + l/2)ln(Kj + 1/2) - (kj - l/2)ln(Kj - 1/2) 



S 



N 

E 

i=i 



(21) 



where > 1/2 (j = 1, . . . , N) are N positive eigenvalues of matrix A", which is a "ratio" of the symmetric 
covariance matrix Q and antisymmetric commutator matrix: 



X = Q^ \ Qjk = ^{qjqk+qkQj), ^jk = {qjqk-Qkqj)- 



(22) 



One can easily verify that if k is an eigenvalue of X , then —k is another eigenvalue. Also, it can be shown 
that all eigenvalues of X are real. It is worth emphasizing that formula ( pl| ) is valid for arbitrary sets of 
operators with c-number commutators (canonical coordinates and momenta, "annihilation" and "creation" 
operators, kinetic momenta and relative coordinates for particles moving in homogeneous magnetic fields, 
etc.). 

In the one-mode case, the eigenvalues of matrix X are equal to ±k (and X^ — k^E2), where 



A = rapp- (xp)^ > = det Q. (23) 

p6l p^.) In this case, different 



(The last inequality is the Schrodinger-Robertson uncertainty relation |^ 
expressions equivalent to formula (|2^) were found in pl| , p8| , p9| . 

Calculating the characteristic polynomial of the 4x4 matrix X in the two-mode case, one arrives at the 
biquadratic equation (for h — 1) [p2| 



I?2K^ + Vn 



0, 



(24) 



where coefficients I?2 and Vq are nothing but quantum universal invariants, i.e., functions which are invariant 
with respect to arbitrary linear canonical (preserving commutation relations) transformations [^0[ : 



P2 = Ai + A2 + 2 {xiX2 P1P2 - P1X2 P2X1) , 



(25) 



= detQ=(p?p2 



~x^ pI {X1P2) 



+2 X1X2 



P1P2' ) \ x1x2~ X1X2 j + {xipi X2P2 - X1P2 x2Piy 

x\ p^ {X2P1) - x^ p^ [xipi) - xf Pi [X2P2) 
2piP2 



PI X1P2 X2P2 + P2 X2P1 XiPi 

-2xTx^ prP2 {xrpi x^2 + xTp2 xWi) 



X2 X1P2 XiPi -h XI X2P1 X2P2 



(26) 
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The symbol A^, obviously, means the combination defined by 
Positive solutions of Eq. (§3) read 



3f) and related to the fcth mode. 



Kl,2 



The reality of ki^2 is ensured by the inequalities 



(27) 



2?2 > 2JVo > 



(28) 



which can be considered as generalized uncertainty relations for two-mode systems (for systematic studies 



of such generalizations see, e.g^ 
Formulae (|2ll), iM) and 



permit us to express the entropic index of correlation (|l|) analytically 
in terms of the covariances of quadrature components for arbitrary Gaussian states. However, the corre- 
sponding expression is very cumbersome, and it is much more complicated than any other entanglement 
measure discussed in the preceding subsections. In the next sections we compare the behaviour of different 
entanglement measures for various concrete physical models. 



3 A three-dimensional cavity with a vibrating wall and two reso- 
nantly coupled modes 

Classical and quantum phenomena in cavities with moving boundaries attracted attention of many re- 
searchers for a long time (see review P2|). Especially popular this topic became in the last decade, being 



known now under the names nonstationary Casimir effect [33| ], dynamical Casimir effect 34 1, or mirror 
(motion) induced radiation One of several theoretical results obtained in the last years was the 

prediction of the exponential growth of the energy of the field under the resonance conditions, when the wall 
performs vibrations at the frequency which is a multiple of the unperturbed field eigenfrequency jl^, ^ . 

A unified description of the field inside an ideal cavity with moving boundaries can be achieved in the 
frameworks of the Hamiltonian approach proposed by Law and developed in (for other references see 
js^ , and for the most recent publications see |39, 4^, Consider a scalar massless field $(r,t), satisfying 



the wave equation ^tt = inside the cavity and the Dirichlet boundary condition $ = on the boundary 
(we assume c = h — 1). We assume that we know the complete orthonormalized set of eigenfunctions (and 
eigenfrequencies) of the Laplace equation V^/Q(r) -I- ^^/^(r) = in the case of stationary cavity. Now 
suppose that a part of the boundary is a plane surface moving according to a prescribed law of motion L{t) 
(for the most recent study of the case when L{t) is a dynamical variable due to the back reaction of the field 
see ^^). Expanding the field $(r, i) over "instantaneous" eigenfunctions fa{r;L{t)), 

<i>(r,i) =^(Z„(0/a(r;L(t)), (29) 

a 

we satisfy automatically the boundary conditions. Then the dynamics of the field is described completely by 
the dynamics of the generalized coordinates qa(t), which, in turn, can be derived from the time- dependent 
Hamiltonian |3^ 

Hit) = ^ E K + ^umhi] + W)T. po^^<^pip (30) 

with antisymmetrical time-independent coefficients 

m^p = -mp^ = lJ dV^^^^^^ffi{v- L). (31) 

For example, in the case of a rectangular three-dimensional cavity with dimensions L^, Ly, Lz, the eigenmodes 
are well known products of sine functions like sinijTkxx/Lx) (or sine and cosine functions in the case of 
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electromagnetic field), labeled by three natural numbers kx,ky,kz, whereas unperturbed eigenfrequencies 
are given by the formula 




-^kM=-\ +^ +^ . (32) 



If one surface of the parallelepiped, perpendicular to the x-axis, moves in the a::-direction (so that the L^- 
dimension of the cavity is a function of time), then p3| 

mkj = (-l)'=^+^^-^4,,„4.,.. (33) 

Jx '^x 

(In the case of electromagnetic field, one should take into account polarizations of the modes, i.e., that fa 
and //3 in Eq. ( ^l|) are vector functions, whose directions are perpendicular, respectively, to the vectors 
{kx,ky,kz) and {ix,jy,jz)- But one can always choose two modes with coinciding polarizations, directed 
along the perpendicular to the plane formed by these two vectors. Then all formulae are the same as in the 
scalar case.) 

We are interested in the case when one of the cavity's walls performs small oscillations with the frequency 
Q, close to the double frequency of some unperturbed mode uj^'^ = 1 (i.e., we normalize all frequencies by 
i^f^), so that the time-dependent frequency uji{i) reads 

cji(i) = 1 + 2ecos(2lZJi), + (34) 

where we assume that ^ 1 and |e| ^ 1. Also we suppose that the unperturbed field frequency spectrum 
includes the frequency ujf'^ = 3 + A with |A| <C 1, but it does not contain frequencies close to bLo'f\ 
A possibility of such a situation was pointed out in ||4^ . An example is a cubic cavity with the pair of 
modes {111} and {511}. Another example is the pair of modes {110} and {510} in the rectangular cavity 
with Lx — V2 Ly (in this case, the common direction of polarization is along the z-axis). Then we have 
two resonantly interacting modes, and it is sufhcient to consider only the part of the total Hamiltonian 
( |3C| ) related to these modes (hereafter we use the symbols Xk, Pj instead of qk, pj for the quadrature 
components of the field, whereas the letter x without indices will mean the usual space coordinate inside the 
cavity) : 

His = ^ {pI +P3) + ^ [1 + 4ecos(2ZUt)] + i [9 + 6A + ecos(2wt)] xj + 3^esin(2cJi) (pixs ~ P3X1) ■ (35) 

The constant parameter /i is proportional to the coefficient mi2 in (|30|). For the rectangular cavity, /i = 
ja;/(12fcj;) if the modes {k^, m, n} and {jx, m, n} are in resonance. Writing ( ^5| ) we have neglected the second 
order terms with respect to e and A. Parameter e has the same order of magnitude as e, but it does not 
affect the solution in the zeroth order approximation . 

Hamiltonian (|3^) results in the following differential equations for the generalized coordinates xi and X3 
(we neglect corrections of the second order): 

ii = - [1 + 4e cos(2wt)] xi + 24fie [cos{2ujt)x3 + sin(2Zjt)i3] , (36) 

X3 = - [9 + 6A + ecos(2Ui)] X3 ~ 24/ie [cos{2ujt)xi + sin(2Ui)ii] . (37) 

These equations have been solved, using the method of slowly varying amplitudes, in [Q. We consider here 
two special cases. 

3.1 Exact (symmetric) resonance 

In the case of exact resonance, 5 = A = 0, the solutions of Eqs. ( |36| ) and ( |3^ ) read 

sin(pr) 



xi{t) = a;i(0) 



Ci cos(pr) + Si 



Pi(0) 



cos(pr) + Cf 



^g^sin(pT) ^^^^^ ^ crp3(0)] , (38) 
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xsit) = 2:3(0) 



TP3(0) 



53+ cos(pt) - C+ 



sin(pr) 



(39) 



Pi{t) = -a;i(0) 



5*+ cos(pt) + C+ 



sin((Or) 



Pi(0) 



C+ cos(pr) + 



sin(pT) 



(40) 



3x3(0) 



53 cos(pt) - C3 



sin(pr) 



+ P3(0) 



-2A^i^2I^ [C3-xi(0)-53-pi(0)], 



C3 cos(pr) - ^3 



_ sin(/9T) 



(41) 



where 



Cj, (t;^) = coshTCOs(fcwt) ± sinhr sin(fcti;f), S^{T;t) = sinhTCOs(fca;i) ± coshr sin(A:a;t), 



1 



^21^- 1, u = 96^2 



(42) 
(43) 



The arguments t ("fast time") and r ("slow time") of the functions C^(r; t) and S^{t; t) can be considered 
as independent variables. Then the following relations hold: 



del 
dt 



dt 



±kC^. 



(44) 



For the modes {111} and {511} of the cubical cavity or {110} and {510} of the rectangular cavity with 
Lx = Ly we have i' = 50/3. Due to this explicit example, we assume that parameter v is large: v ^ 1. 

Symbols Xk and pk in equations (|3|)-(41) can be considered both as classical variables and quantum 
operators in the Heisenberg picture, due to the linearity of the problem (or due to the quadratic nature 
of Hamiltonian (|30|)). Using equations (^-(^, one can calculate mean values of squares and products of 
canonical variables (operators) at any moment of time, provided such mean values were known at the initial 
moment t = 0. We confine ourselves to the case when initially the field modes were in thermal states with 
the mean photon numbers (f?i — l)/2 and (^3 — l)/2, where 9k = coth(fc/?/j/2), f3k being inverse absolute 
temperature in dimensionless units. In the natural case of equal initial temperatures of the modes, the 
following relations hold: 



^13 



1 > 



'31 



1 

> -. 

- 3 



(45) 



The normalized mean energies in each mode, £k — -r ^k-^k 

the definitions of the covariance entanglement coefficients dl8) and (f9|)), depend on time as follows |4 



{p1 + ujIxI) / {2ujk) (namely these quantities are used in 



£1 



cosh(2T) 



sin^(pT) 



{1 + 2v 03i) +cos^ (pr) 



- sinh(2r) 



dn(2pr) 1 



(46) 



£3 = Y] cosh(2T) 



sin^(/9T) 



(1 + 2i/6'i3)+cos2(pt) 



- sinh(2T) 



sin(2pT) I 



P 



(47) 



Calculating the covariance entanglement coefficient, one should use, instead of variables Xk and pk, the 
normalized variables Xk — y^uJkXk and pk — Pkl (in our case uj^ = k): see equation (po[). After some 
algebra we have obtained the following expressions: 



3^ = 



4£l£3' 



y = 



^1+^3 



(48) 
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T = sin (pr) < cosh(4T) 



The determinants of the covariance matrices for each mode have been calculated in 



(49) 



detQn = -^Ofgl 



gl = cos4(pr) + sin2(2pr)^^|^^ + sin4(pr) 



2v( 



For the first mode, 

2 



'31 



1 



2v- 1 



(50) 



For another excited mode one should interchange indices 1 and 3 in (jsy). Since the evolution of the total 
system is unitary in the case discussed, the total determinant does not depend on time: det Q = O^O^/IQ. 
Therefore the purity entanglement coefficient has the form 

>C = l-(5iff3)"'. (51) 
Eqs. (^fj) and ( p3| ) lead to the following explicit formula for the entropic entanglement measure Ic (|l|): 

Ic^\Y. + 1) In (0,3z + 1) - (^?.5* - 1) In - 1) - + 1) In (^, + 1) + (^.-1) In (&,-!)]. (52) 

i=l,3 

We see that despite the exponential (although non- monotonous in the high- temperature case 6*^ ^ 1 [Q) 
growth of energy of each mode, all entanglement coefficients exhibit strong (quasi)periodic oscillations as 
functions of the "slow time" r, going to zero when pr = nn. 

In the simplest case of the initial vacuum states of each mode {6i — 9^ — 1) we have 



4i/ 



(2i/- 1)2 



2 2 2 

5i = 53 = % 



sin'*(pr) cosh(4T), 



8i/ 



8;^ sin^(/9T) 



(2:/- 1)2 + 8j^sin''(pT) 



/, = (5o + l)ln(.go + l)- (.9o-l)ln(.go-l)-21n2: 



;sin''(pT), 
sin'*(pr) 



In 



sin {pr) 



The approximate equalities in ( p35D and (|5g) hold for ly ^ 1. Under this condition, £i w £3 
that for r > 1 we obtain 

37- ./Isin2(pr) « VZ. 



(53) 

(54) 
(55) 
(56) 

i cosh(2T), so 



The evolution of functions C{t) and 3^^(t) for the initial vacuum case is shown in Figure 

For high-temperature initial states (6*1 ^3 ^ 1), the entanglement coefficients do not depend on the param- 
eter 1/ (if v ^ 1) for almost all instants of time, more precisely, under the condition | cos(pr)| ^ p~^ ~ i/~^/2; 



C 



sin2(2pr) (031+013-2) 
4 + sin2(2pr) (^31 + ^13 - 2) 



Ic = In (3153) ~ In 



l + -sin2(2pT) (031+^13-2) 



(57) 



(58) 



(The last approximate equality in (^7|) holds for r > 1. The simple formula for Ic is obtained in the limit 
case 01.3 00; there are some corrections of the order of 0f 3 for finite initial mean numbers of photons.) 
For the maximal possible value of the coefficient 0i3 = 3 (in the case of true initial thermal equilibrium), 
the maximum values of the expressions (|57|) and (|58|) (which are achieved when sin^(2pr) — 1) are equal 
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to Umax = 1/4 and j^™"^) = lii(4/3) « 1/3. Note that in the hmit high-temperature case, the purity 
entanglement coefficient coincides identically with one of possible forms of the "compact entropy" (another 
compact parameter, tanh(/c), was introduced in p^) 



Jc = 1 - cxp (-/c) . 



(59) 



Figure |^ shows the evolution of entropic entanglement measure IdT) for vacuum and high-temperature 
initial states. Note that one has Oi ~ 140, if Lq — 1cm and T = 300 K. For this value of 6*1, the plot of the 
compact entropy JAj) becomes indistinguishable from the plot of the purity entanglement coefficient Cij). 
The functions £(t) and 3^^(t) are compared in Figure ||. We see that two functions are very close in some 
intervals, although their maxima are different (because the value p ~ 5.7 is not very large for the chosen 
parameter v ~ 50/3). 

In the high-temperature case, intermediate nonzero minima of the entanglement coefficients (besides 
exact zero minima at the instants r„ — njd p) are observed at the moments of "slow time" when the modes 
approximately exchange their purities ||44| . The positions of these additional minima for C and Ic are 
determined by the condition cos(/3t) — 0, so that 



r 



For ^ 1 we have 



(6>3i -t- 6>i3 + 2) 
4r/2 + 1 + 2v {031 + 613) ' 



& 



j(rnin) 



731 



= In 



713 



1 



2z/| 



731 



713 



■2) 



(2^.-1)2 



2z/ 



(60) 



(61) 



On the other hand, the intermediate minima of y are much smaller. Indeed, the minimum of the expression 
inside figure brackets in Eq. (MS) is achieved for (neglecting corrections of the order of p~^) 



tan 2pT = -tanh(4r)-i 

p &1-&3 



At this moment of time we obtain 



Av cosh(4T) 



4£i£3 w 6163 cosh2(2T), 



so that for r > 1, 



-4r 



731 



X w 2e 



-4t 



and it is clear that the intermediate minimum of y does not exceed the value 3^*. 

Therefore, we arrive at rather paradoxical situation, especially for realistic values of parameters u and 
9i^3. According to Figure |[ the intermediate mimimum value of ^-coefficient in the high-temperature case 
is only twice less than the maximal value. Moreover, this high-temperature intermediate mimimum value 
is bigger than the maximum value in the vacuum case (see Figure |l). Thus, the /^-coefficient tells us that 
for cos(pr) = 0, two modes are "more entangled" in the high-temperature case than in the case of initial 
vacuum state (or at least have the same order of entanglement, according to the /c-coefficient in Figure 
whereas the covariance entanglement coefficient y shows that two modes become practically disentangled at 
this instant of time. 

The resolution of this "paradox" is as follows. According to Eqs. (^8|), (^) and (^), the function T gives 
the upper limit for squares of any elements of the "off-diagonal" block Q12 of the covariance matrix Q (||), 
whereas functions £k give the bounds for the elements of "diagonal" blocks Qkk- This happens because T 
and are based on traces of the covariance submatrices. Therefore, if 3^ ^ 0, this means that all elements 
of matrix Q12 responsible for the intermode correlations (at least for the Gaussian states considered in this 
section) become negligible in comparison with the variances x^Xk and PkPk of the quadrature components. 
From the physical point of view, it is equivalent to disappearance of correlations between the two subsystems, 
i.e., their disentanglement. 



10 



On the other hand, the coefficients L and Ic are based on determinants of the covariance submatrices. 
But it is well known that the determinant of a matrix can be quite small even if all elements of the matrix are 
big, and this is the reason of the qualitative difference in the behaviour of the "covariance" and "entropic" 
entanglement coefficients. This is clearly seen from the last expression in Eq. (H), which shows that the 
value of the purity entanglement coefficient C depends of the matrix R — Q12Q22 Q'^iQii ■ Using easily 
verified formula det(£' + a) ~ Tra, which holds provided all elements of matrix a are small with respect to 
unity, we can simplify formula ^ in the case of small entanglement as follows: 

C^Itt{QuQ22Q2iQii)- (62) 

But each matrix Q^^ (fc = 1, 3) contains the denominator det Qkk, which can be much less than any element 
of matrix Qkk- If this happens, then the inequality Tr {Qi2Q2i Q21Q11) > Tr (Q12Q21) / (TrQiiTrQ22) 
becomes quite possible. Just such a situation takes place in the example considered. Although diagonal 
elements XkXk and PkPk of matrices Qkk grow exponentially with time, these matrices have also exponen- 
tially growing off-diagonal covariance elements XkPk (this means that each mode occurs in highly- correlated 
quantum state p6t with quadrature correlation coefficient r = XkPk/ {xkXk PkPk)^^'^ approaching the unit 
value), so that det Qkk does not grow unlimitedly with time, exhibiting only relatively small oscillations. For 
this reason, elements of matrices Q^^ have the same order of magnitude exp(2T) ) as elements of matrices 
Qkk themselves. On the other hand, elements of matrix Q12 have an order of exp(— 2r) at the moments of 
intermediate minima. Therefore, the exponential time dependences are canceled in the measures based on 
determinants, resulting in the inequalities j7c,>C ^ for cos(pt) ~ 0. 

This example permits us to make a conjecture that the covariance entanglement coefficient y is not only 
simpler from the point of view of calculations, but it could be preferable from the physical point of view, 
because it is more sensitive to entanglement than entropic and purity measures. Other arguments in favour 
of y can be found in [jis) . 



3.2 Asymmetric resonance 



An interesting feature of the Hamiltonian (|3^) discovered in |Q is a possibility to compensate one detuning 
(e.g., S) at the expense of another. In particular, an exponential growth of the energies of both modes can 
be obtained under the conditions of "asymmetric resonance" 



5 = e, 36- €1^/2. 
In this case the quadrature components depend on time as follows: 



(63) 



Xi 



it) = 



.Tl(0) 
, ^3(0) 



1 - - Cf {2Rt; t) + - cos 



[Cf (2i?r;t) -cos(/)i] 



P3(0) 
12/i 



[S^{2RT;t) + sm(j)i] , 



1 - - S^{2RT\t) - - sin 01 



(64) 



xs{t) 



^3(0) [I 

, ^i(O) 
12/i 



1 - ) cos 03 



-C^{2RT;t) 



1 I sin 03 



--S^{2RT-t) 



[C^{2RT-t) 



- cos 03 - 



pM rc- 
l2^l 



[S^{2RT-t)- 



' sm 03 



(65) 



Pi{t) 



-xi(0) 

X3(0) 

" 4/x 



1 - - I St(2RT;t) + - sin 01 
J/ / J/ 



Pi(0) 



2\ 2 
1 - - Ct{2RT;t) + - cos 01 

V I V 



[S'i+(2i?r;<) -sin0i] 



P3(0) 
12^ 



[Ci+(2i?T;i) -cos0i] 



(66) 
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xijO) 



1 - - ) sin03 + -St(2RT;t) 



2\ 2 
1 - - cos(/)3 + -Ct(2RT;t) 



where 



[53+(2i?r;i)-sin(/.3] 



ipkiT', t) = kiot — 2Jr, 



p_m 



[C+i2RT;t)-cosh] , 



(67) 



R= 1 



and all terms of the order of ©(i^^^) have been neglected (as well as the corrections of the order of 5 ^ e in 
the amplitude coefficients). 

The (normalized) mean energies of each mode depend on time as follows E^]: 



4\ 4 
1 - - j cosh(4i?T) + -iP{t) 



— [cosh(4i?T) + 1 - 2V' (r)] , 



2 



4 4 
1 - - + -V^r) 



+ — [cosh(4i?T) + 1 - 27/- (7 



(68) 
(69) 



where 



i/'(t) EE cosh(2i?T) cos(2 Jt). 



The energy of the third mode is significantly less than the energy of the first mode, if 3> 1. For this reason 
this regime of excitation was named "asymmetrical". For r > 1, £3 /Si « Q/v. Note, however, that for the 
cubical cavity with v = 50/3, the energy of the third mode is only three times less than that of the first one. 
It is important, nonetheless, that the rates of increase of the energies of each mode are almost twice bigger 
than they were in the case of the strict resonance discussed in the preceding subsection. 

The covariance entanglement coefficients can be written again in the form (48), but with £1,3 given by 



( p^ ) and (|6^) . The function T in the asymmetric case reads (neglecting corrections of the order of ^ with 
respect to the main terms) 



T ^ 2 



cosh^(4i?T) + sinh^(2i?T) - cosh(6i?T) cos(/)o 



+2i/"i [cos0o{cosh(2i?r) +3cosh(6i?T)} - 2 cosh^ (2i?T) cos^ 0o - 2cosh(4i?T) - 2 sinh2(4i?r)] 



-6I3 cosh^(2i?r) - cosh(2i?T) cos 0o 

-2v-^ [cos(/)o{cosh(6iT!T) +3cosh(2i?T)} - 2 cosh^ (2i?T) cos^ 0o - 2 cosh(4i?T)] ) 
-26I16I3 ( cosh(4i?T) [cosh(2i?T) cos (/)o - 1] 



-2v ^ [— 2cos(/)o {cosh(6i?T) + cosh(2i?T)} + 2 cosh^(2i?r) COS^ 00 + 

4cosh''(2i?T) - 2] 



(70) 



where — —2 Jr. 

If r — > 00, then (for v ^ X) 



fl2 

2v ' 



£1 ~ ^ e , 



£3 « l^e^^^ 
2v 



so that 3^ ^ 1, whereas y Consequently, the coefficient y is preferable when the energies of 

subsystems are essentially different. 



The purity and entropic entanglement coefficients are given by Eqs. ( |51| ) and (52), with 

glir) = ! + -[(!- ^31) V^r) - I + 631 cosh2(2i?r)] (71) 
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and (73 obtained from (|7l| ) by means of the replacement 1 <-> 3. We see a significant difference from the strict 
resonance case: now functions 5i,3(t) increase exponentially with time for r 3> 1. Asymptotically, each 
mode appears in a highly mixed quantum state, with det Qn = det Q33 = 9i93 exp(4i?r)/(2z^). The purity 
entanglement coefficient ( |5l| ) tends asymptotically to the unit value independently of the initial temperature 
(or coefRcients 9k)- 

£ w 1 - ^exp(-4i?r), r > 1. 

For r 3> 1 the entropic entanglement coefficient grows unlimitedly: 1^ ^ lii((7i53) ^ 4_Rt. Therefore in 
Fig. H we compare the compact parameter J'c{t) ( p9| ) with the functions 3^(t) ([l^), £ (^, and [£(t)]^/^ 
for the initial vacuum state. Since all formulae in the asymmetric case are obtained neglecting terms of the 
order of i^^^, we use the value v — 100 in the illustrations. The difference between the asymptotical values 
of the functions for t ^ 1 and the correct value 1 shows the accuracy of approximation (about 2%). The 
dependences of the entanglement covariance and purity coefficients y and £ on the "slow time" r for the 
initial vacuum and high-temperature state are shown in Fig. ^ Remember that in the high-temperature 
case the coefficient £ tends to the compact entropic coefficient jTc- 



4 Fabry-Perot cavity with an oscillating boundary 



The problem of the scalar massless field in a ID cavity formed by two infinite ideal plates whose positions 
are given by x/e/t = and 



'bright 



L{t) = Lq{1 + s sin [pujit]) , |e| < 1, Wi = ttc/Lq, p = l,2,. 



(72) 



was solved in |Q. The only component of the operator vector potential of the electromagnetic field A{x,t) 
in the Heisenberg representation can be written as 



A{x,t) = \bn^<-^Hx,t) + h.C. 



bn, bl 



— 5nk 



(73) 



where 





irkx 


sin 


[Lit] 


1 






2 





{pi")(r)e--*-pL1(r)e^-*} 



UJn ncji. 



(74) 
(75) 



The normalization factors 2/y/n in ( |7^ ) are chosen in such a way that the energy of the field in the stationary 
case can be represented as a sum of energies of independent mode oscillators. The coefficients p^^\t) satisfy 
an infinite system of coupled equations (fc = ±1, ±2, . . .; n = 1, 2, . . .) 



d7^^ 



(") _ 



(n) 



(76) 



which was solved in |Q (here we confine ourselves to the simplest special case of solutions found in 
corresponding to the strict resonance). 



Due to the initial conditions p^,"'' (0) — Skn the solutions to ( [76[ ) satisfy the relation p 



The non-zero coefficients pin^ read JlS] 

Pj+mp' (''') 



(fe+np) 
j+mp 



where 



r{l + n + j/p) (ctk)"-™ 
r (1 -|- m + j /p) r (1 -|- n — m) 

K = tanh(pT) 



F (^n + j Ip , — m — j /p ; 1 + n — m ; k^) , 



if k. 

(77) 
(78) 
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and -F(a, 6; c; z) is the Gauss hypergeometric function. The functions (|77| ) are exact solutions to the set of 
equations (|76| ) relating the coefficients with different lower indices. Besides, these functions satisfy another 
set of equations, which can be treated as recurrence relations with respect to the upper indices ||45|| 



dr 



o*") = n 



Jn-p) (n+p) 



d 



(p-n)* _ (p+n) 
r—ni rm 



}, n>p, pW^o 

} 



n = 1,2, 



The consequences of equations ([Zq), ([79|) and (|80|) are the identities 

oo 

'^P'm''*Pm = nSnk , n,k = 1,2,. 



m— — oo 



OO 



) (n)* (n) 



n=l 



oo _j 



= (5„y , TO,j = l,2,. 



0, m,j = 1,2,... 



(79) 
(80) 

(81) 
(82) 
(83) 



We suppose that after some interval of time T the wall comes back to its initial position Lq. For t > T, 
the field operator assumes the form 



2 

x,t) —} sin (nnx / Lq) \a„ 

n—l ^ 



(84) 



where operators am are related to the initial operators bn and bl^ by means of the Bogoliubov transformation 
(tt = ^SLJiT) 

OC j 

Im 



n=l 



m = 1,2, . 



(85) 



The commutation relations |^d„ , a|,j = Snk hold due to the identities (^l|)-(p3|) which are nothing but the 
unitarity conditions of the transformation (pq). These commutation relations together with the expression 



for the energy of the field 



H 



8^ Jo 



dx 



dA 
'dt 



dA 
dx 



E '^^ ( 



n=l 



(86) 



convince us that a„ and a|j are true photon annihilation and creation operators a.t t > T (like the operators 
bn and bl^ were 'physical' ones at t < 0). 

4.1 Intermode entanglement in the parametric resonance case {p = 2) 

Our first goal is to calculate the entanglement coefficients between different modes in the case of the para- 
metric resonance, p = 2. If the initial state of the field was vacuum with respect to the initial operators bn'- 
bn\0) = (we use here the Heisenberg picture), then the covariance entanglement coefficient between the 
rth and sth modes is 

1 1/2 

|(a^a,)|2 + \{alas)\^ 



yr.s = 



(^{alar) + 1/2) ((ala.) + 1/^ 



(87) 
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Using jssl), one can express the average values contained in formula js^ ) as (assuming hereafter lui = 1) 

oo ^ oo ^ 

(a,a,) = ~V^sJ2 = -V^E -pi'^Vt^*, (88) 



(alas) = vrs 



oo 



oo _j 

•E- 



(89) 



where the coefficients should be taken at the moment T, thus their argument is tt- Strictly speaking, 
the expressions in (^) and (|8^) have physical meanings at those moments of time T when the wall returns to 
its initial position, i.e. for T = Ni: /p with an integer N . Consequently, the argument tt of the coefficients 

pf^m ™ (H) (H) assumes discrete values t*^^) = NeTr/{2p). One should remember, however, that 
something interesting in our problem happens for the values r ~ 1 (or larger). Then ~ 3> 1, and the 
minimal increment At ~ e is so small that tt can be considered as a continuous variable (under the realistic 
conditions, e < 10^^ [l^ ). For this reason, we omit hereafter the subscript T, writing simply r instead of 
Tt or T^-'^K 

Differentiating the right-hand sides of equations ( ^8|) and (^ ) with respect to the 'slow time' r, one 
can remove the fraction 1/n with the aid of the recurrence relations ( [tq]) and (^0|). After that, changing if 
necessary the summation index n to n ± p, one can verify that almost all terms in the right-hand sides are 
cancelled, and the infinite series are reduced to the finite sums. For p = 2 we obtain the equations (taking 

J") 



(90) 



into account that all functions pin are real in the strict resonance case, according to Eq. {u7h ) 



d(a,.as)/dT 



d{ala,)/dT = Vfi U'V-] + , d{alar)/dT = 2r p^^V-r 



(91) 



For p — 2, only odd modes can be excited from the initial vacuum state. In this case, the hypergeometric 
functions in the formula (|77| ) for coefficients pr"'' with j — I are reduced to some combinations of the complete 
elliptic integrals of the first and the second kinds 



tt/2 



da 



sm a 



,11 

= -F - , - : 1 ; K 



E( 



r^"^ I TT / 1 1 

k) = y daV 1 - sin^ " 2 ^ ( ^ 2 ' 2 ' ' 



so that equations (^) and ( ^ ) can be integrated for any values of r and s: see 46 or Appendix ^ for 
technical details. In particular, for the first few modes we find 



{k^Y? - 2EK + E^] , 



[k2(4 - k2)k2 - 2(2k'' - 4)EK + (4k'* - 



■4)E2 



2^/3 



2^3 



[^2k2 - 2EK + (1 + e2] 



2)EK + E" 



2V5 



[^2(^2 + 8)K2 - 2(k^ + 8)EK -h (8k^ Tk^ 



— (2k2 + 1)K2 + -(2^4 - 2^2 - 3)EK + (k^ 
5 5 



)E^ 



1)E2 



(92) 
(93) 
(94) 
(95) 

(96) 
(97) 
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(agas) = [k2(k2 + 2)(/«2_ 2)K2 + 2(2^6 2^2 + 4)EK 

-4(k2 + i)(k4-k2^1)E2] , (98) 



(a. 



3°5 ) = -^S^ [^^{Tn^ - 4)K2 + 2(8/c4 - 15^2 + 4)EK - {Ak^ - Wk^ + 4)E2] , (99) 

2 



457r2K 



£1 = — K 2E-K^K , (100) 
= [(3^2 - 2) K (2E - K^K) + 2 (l + k^) e^] (101) 

£5 = - ^, . [k^(47k^ - 30k^ - 8)K^ + 2(4k^ - 47k^ + 26k^ + 8)EK 
457r'^K* 

-2(^2 + 1)(4k'' + llK^ + 4)E2] , (102) 



where k = -\/l — a^d we used £r = (ala,.) + 1/2. 

In Figure |^ we show 3^i^3 and 3^3,5. We see that the entanglement is strongest for the lowest modes. 
However, for any pair r, s the coefficient 3^r,s tends asymptotically to the unit value when k — > 1. To prove 

(n) 

this property, one should use the asymptotical forms of the coefficients pm for r — > 00, i.e., for k 1. 
Namely, replacing the hypergeometric functions in (|7^) by their values for the unit argument |47[ , 

u u . r(a + 6+l) 

^("'^'" + ^ + ^'^) = r(a + i)r(6 + i) ' 

one obtains the following asymptotical formulae (see also [M): 



P^:,UrM=pLV,(r) = ^^i^, r^oo. (103) 



Consequently, for r 1 we have 

(ala,) « -{ara,) « ^^(-l)(--^)/2 + 0(1), (104) 



and the leading terms in the numerator and denominator of the fraction in ( |87| ) become the same for r ^ 00. 

In the case of detuning from the strict resonance, characterized by some dimensionless detuning parameter 
7, the coefficients pm"* become complex. However, their asymptotical forms differ from (103) only by some 
phase factors . Since the covariance entanglement coefficient (^) depends on the absolute values of the 
second-order moments (ajd^) and (0,^03), these phase factors do not influence the final result, namely, that 
y^s — s- 1 when T —^ 00, unless the dimensionless detuning parameter exceeds the critical value 7 = 1, when 
the generation of photons from vacuum becomes impossible. 

In the case of the initial vacuum state, the state of the field at the subsequent moments of time remains 
Gaussian p^ , and the purity entanglement coefficient can be calculated by means of formula (|8|). In a 
generic case, the determinant of the symmetrical 4x4 matrix Q (|2^) contains 17 different terms. However, 
in the specific case involved all covariances between the "coordinate" and "momenta" operators turn out to 
be equal to zero identically: XkPj = 0, and for this reason the determinant of the covariance matrix for the 
ith and jth modes can be factorized in the following simple form: 

detQ = {(^p.p^^PiPj - ^l,p^ (yO-x,x,crxjXj - f^l^x^ ■ (105) 
Nonzero covariances are given by the following expressions: 

cTxiXj = XiXj = ^{ajaj + a]dj) + Re(didj), ap^p. = ^ ^("loj + ajoi) - Re(aidj). (106) 
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Introducing the correlation coefficients, 



' Pipj 



\/^PiPi^': 



PiPi 



we can represent the C (H) and Z (|T^) entanglement coefficients between the ith. and jth modes as 

r — 1 - 



1 - r2 



PiPj 



1 - r2 



PiPj 



\{^-\rl.x^){^-\rQ 
If all correlation coefficients are small (in particular, if t ^ 1), then 



1 - r2 

P^P] 



(107) 



(108) 



(109) 



PiPj 



When r — > oo, then, due to equations (104), (10£) and (107), the coefficients (Jp^p^ linearly grow with time 
in such a way that the momentum correlation coefficient Tp.p. tends to the unit value. At the same time, the 

coefficients UxiXj and Tx^xj tend to some finite limit values. Therefore, the purity entanglement coefficient C 
and the distance entanglement coefficient Z approach the unit value. Using the asymptotical formulae for 
the complete elliptic integrals E^, 



4 1/4 
K(k) « In - + - In - 

K 4 V K 



1 + 



E( 



i+ifml 

2 I k 



one can see that for r 1, 1 — £ ~ 1 — Z ~ r '^Z'^. In particular 

0.88 



1-C 



13 



44 
577 



V 3r 



^219 



1.19 



Calculating the entropic entanglement measure (|i|) one should take into account that the reduced entropy 
of any two-mode subsystem depends on time in the case involved (in contradistinction to the case considered 
in the preceding section), because the evolution of each finite-dimensional subsystem is not unitary. This 
entropy is determined by two eigenvalues of the corresponding 4x4 matrix QQ,~^, which are given by 
formula (p7|). The reduced entropy of the fcth mode is determined by the single number 

fk = y/<^PkPk'^XkXk , (110) 



as soon as the coordinate-momentum covariances are equal to zero in the strict resonance case considered. 
The explicit formula for the entropic entanglement measure between the fcth and nth modes becomes (for 
the initial vacuum state of the field) 



^ [ (/, + 1/2) In (/, + 1/2) - (/, - 1/2) In (/, - 1/2) 



j—k,n 



[ {fL + 1/2) In {fL + 1/2) - {fL - 1/2) In {fL - 1/2) 



5=±1 



where 



2/L = 



PkPk XkXk + PnPn XnXn + 2pkPn XkXn + 2 {pkPk PnPn " PkPn'^) {xkXk 

p^ x^ + S;;^;: + '2pkP^ XkX^ - 2 {wPk P^ - PkP^^ 



1/2 



(111) 



1/2 



(112) 



The behaviour of different entanglement coefficients is shown in Fig. 0. All of them monotonously tend to 
unity with the course of time, but much more slowly than in the case of asymmetric resonance in the 3D 
cavity (due to interaction with other resonant modes). 
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4.2 Entanglement in the "semi-resonance case" (p = 1) 

A qualitatively different behaviour of all characteristics of the field is observed in the "semi-resonance case" , 
when the frequency of the oscillations of the boundary coincides with the fundamental field eigenfrequency 
(p = 1) 48 . In this case one should put j = in formula ([77|), and all coefficients pm^ with negative 



lower indices m are equal to zero identically. As a consequence, no photons can be created from the initial 
vacuum state, which is clearly seen from equation (p9|). If initially the field was in non- vacuum state (at least 
for some mode), then the total number of photons in all modes is conserved, although the total energy grows 
exponentially due to "heating" the high-frequency modes (at the expense of "cooling" the low-frequency 
modes). 

We suppose that initially only the first mode was excited, while all the others were in the vacuum state. 
Then the dynamics of all modes is described by means of the unique coefficient 

= (tanhT)™-VcoshV. 

If initially the excited mode was in a coherent state, then all second-order central moments connecting 
different modes are equal to zero, resulting in zero covariance entanglement coefficient: 3^^°'' = 0. However, 
for other initial states we obtain nonzero values of 3^r-,s- 
If initially the first mode was in the Fock state \n), then 

y^r'^ = ''^''^^ , (113) 

v/2« + l/2) « + i/2)' 



where 



^ (^) n- (tanhr)™-! ^ 



cosh^ 



T 

t2 t2\ 



If initially the first mode was in a squeezed vacuum state — exp i?(6| — 6f )/2 |0) with the average 
number of photons = sinh^(i?), then 



ysqz ^ CrCsVt^l(2t^l + l) 



If initially the first mode was in an even/odd coherent state 

|ai) ± I - ai) 



|a± > = 



V2[l±exp(-2|ai|2)] 
then the mean numbers of photons are given by the formulae 

= laiptanhdaip), = cothdaip). 

In both cases, the entanglement coefficient can be written as 



For big enough number of photons in the initial squeezed and even/odd states, i^i >• 1, the entanglement 
coefficient becomes very close to the maximal possible unit value, ii ^ ^ 1, but with increase of time y 
eventually goes to zero, because Cr.sif) ^ for t ^ oo. In the case of the initial Fock state, the maximal 
value of y does not exceed 1 / y/2. A typical behaviour of the covariance entanglement coefficient between the 
first and second modes for the initial Fock and squeezed states is shown in Fig. ^. The behaviour of 3^m,n 
for the initial thermal and even/odd states is very similar, especially for large mean numbers of photons. 
The evolution of the mean number of photons in the first and second modes is shown in Fig. |9[ 

The momentum-coordinate covariances turn out to be equal to zero again (as in the case of p = 2), 



therefore we need only two correlation coefficients defined in (107), in order to calculate the purity and 
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distance entanglement coefBcients (in the case of the initial squeezed state of the first mode) with the aid of 

(117) 



Eqs. (108) and (|109|). These correlation coefficients are as follows, 

^ _ xC.(r)0(r) ^ _ AC. (r)O(r) 



[1 + xQl{r)] [1 + xC|(r)] 01-ACf(r)] [l - AC|(r)] 



where 



X = e2«-1, A = l-e-2«. 



The time dependences of the C and Z entanglement coefficients are compared in Fig. We see that the 
full and dashed curves are very close, especially for large mean numbers of photons. 

5 Conclusion 

We have compared time dependences of several functions characterizing the degree of entanglement between 
field modes of ideal cavities with resonantly vibrating walls for different models of such cavities. All these 
functions (the "standard" entropic entanglement measure for Gaussian states, covariance entanglement co- 
efficient introduced in |l^, |l5| , distance entanglement coefficient introduced in p3[ | , and purity entanglement 
coefficient introduced here) are based on the second-order covariance matrix of the field quadrature compo- 
nents. In spite of having different analytical forms, the coefficients concerned show similar qualitative (and 
in certain cases even quantitative) behaviour for each fixed model. Therefore, the covariance entanglement 
coefficient, being the simplest from the point of view of calculations, seems to be the most convenient, es- 
pecially compared with the entropic entanglement measure, whose calculation requires tremendous efforts, 
giving practically the same information on the degree of entanglement. Moreover, an example at the end of 



section 3.1 shows that the covariance entanglement coefficient (based on traces of covariance submatrices) 
can be more sensitive to entanglement than other measures (which are based on determinants of covariance 
submatrices) . 

On the other hand, the behaviour of each selected entanglement coefficient turns out to be completely 
different for different kinds of cavities. For the three-dimensional cavities with accidental degeneracy of the 
spectrum of eigenfrequencies, the entanglement coefficients exhibit oscillations in the case of "symmetric" 
resonance, remaining relatively small for all instants of time. Moreover, they go to zero periodically, despite 
that the energy of each mode increases unlimitedly. In the case of "asymmetric" resonance, fast (in the 
"slow time" scale) oscillations of the entanglement coefiicients are also observed, but all these coefficients 
tend to the maximal possible unit value with increase of time. For the model of one-dimensional (Fabry- 
Perot) cavity with equidistant spectrum, all entanglement coefficients monotonously go to the unit value 
in the parametric resonance case. In the "semiresonance" case, they rapidly reach the values close to 
unity and remain at this level for some interval of time (which increases with increase of the initial mean 
number of quanta), but eventually they decay to zero. Therefore, this study adds some new features to our 
understanding of the behaviour of fields in cavities with vibrating boundaries, in addition to results obtained 
earlier in |l|, g, |45[|6[ |0) . 
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A The Bogoliubov coefiicients in the ID parametric resonance 
case 

The nonzero coefficients P2m+i '^^ parametric resonance case {p — 2) read Q 

fi1 (-l)"T(m+ 1/2)k™ / , , 9x ^ 

= r(l/2)r(l-fm) -^(^ + V2,-l/2;l + m;..^), (A.l) 



19 



(1) (-l)"T(m + l/2)r(3/2)K"+i , , , ,x , 

= ^r(2 + m) F(m+l/2,l/2;2 + m;«2)^ (^2) 



In particular {R, = Vl — k^), 

p['^ = -E{k), p^1\ = — [^{k) - k^n^^)] , (A.3) 

4" = [(1 - 2^') E(«) - S:^KM] p«=-3^[(2-^2)E(^)-2^2K(^)], (A.4) 
= _H_ [(8^4 _ 3^2 _ 2) E{k) + (-4^4 + 2k^ + 2)K(k)] , 
p^^l = - — ^ [(2^4 + _ 8) E(k) - (k^ + _ 8)K(k)] . 
The general structure of the coefficients P2m+i terms of the complete elliptic integrals is GgI 



P2I+1 = — ^ [/™ E(^) + K(k)] (A.5) 



P-'L-i = -^^i hn E(k) + les^ {k") Kin)] (A.6) 

where fmix)^ gmix),Tm{x), Smix) are polynomials of the degree m which can be found from the recurrence 
relations (uw. 



B Calculation of integrals 



To calculate, for instance, the average value {a^a^), we use equations (pl|), (A.3) and (A.4), replacing the 
derivative over r by the derivative with respect to k in accordance with the relation dn — 2k^dT. In this 
way we arrive at the equation 

= [(1 + E'('*) - 'i'K2(/^) - 2^'~k'E{k)K{^)] . (A.7) 

Taking into account the differentiation rules 

dK(K) E(k) K(k) dE(K) E(k.)-K(k) 



dK KK^ K ' dn 



(A.; 



we may suppose that the factor in the denominator of the right-hand side of equation (A.7) comes from 
the derivative dK/dK. Thus it is natural to look for the solution in the form 

{a\a^ ) = 1^ [A(«:)K2(k) + B(a«)K(k)E(k) + C(k)E2(k)] , (A.9) 



where A{k), B{k) and C{k) are some polynomials of k. Putting the expression (A.9) into equation (A.7) we 
obtain a set of coupled equations for the coefficients of these polynomials, which can be resolved recursively. 
The equations for other second-order moments can be integrated in the same manner. 
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Figure 1: The covariance entanglement coefficient squared (thick hne) and the purity entanglement 
coefficient £ (thin line) versus "slow time" r for two interacting modes {1, 1, 1} and {5, 1, 1} in a 3D cubical 
cavity — 50/3) under the condition of strict ("symmetric") resonance and for the initial vacuum state. 

Figure 2: The entropic entanglement measure Ic ( |5^ ) versus "slow time" t for two interacting modes {1,1,1} 
and {5, 1, 1} in a 3D cubical cavity {v = 50/3) under the condition of strict ("symmetric") resonance, for 
the initial vacuum state (thick line; = ^?3 = 1) and high-temperature state (thin line; 9i = 3^3). 



Figure 3: The functions C{t) (thin line) and 3^^(t) (thick line) for two interacting modes {1,1,1} and 
{5, 1, 1} in a 3D cubical cavity (z^ = 50/3) under the condition of strict ("symmetric") resonance and for the 
high-temperature initial state with 9i = 140. 

Figure 4: Time dependences of different entanglement measures under the condition of "asymmetric reso- 
nance" (S3h, for ly =^ 100 and the initial vacuum state. Thick line: the covariance entanglement coefficient 



3^(t) (lis). Thin lines from bottom to top: the purity entanglement coefficient £ (51), the compact entropic 
entanglement measure JdT) (M), the function [£(t)]^/^. 



Figure 5: The covariance entanglement coefficient 3^(7^ (|T^) and the purity entanglement coefficient C{t) 
( |5l| ) under the condition of "asymmetric resonance" (^) with i/ — 100, for the initial vacuum state with 
61 = 9^ = 1 (monotonous dependences) and high-temperature state with 9i = (oscillating functions). In 
both cases, upper curves correspond to 3^(t) and lower curves correspond to C{t). 

Figure 6: The covariance entanglement coefficient Xi,m (^T]) ™ resonance (p = 2) cavity versus the 

compact parameter k = tanh(2r) for the vacuum initial state. Full curve: 3^1,3; dashed curve: 3^3,5. 



Figure 7: Different coefficients characterizing entanglement between the first and third modes in the ID 
resonance {p = 2) cavity versus "slow time" t (in the insertion) and the compact parameter k = tanh(2T), 
for the vacuum initial state. The order of the curves in the main figure, from top to bottom: covariance 
entanglement coefficient y (^); compact entropic coefficient Jc (|5^); purity entanglement coefficient C 
(|5l|); the square of the covariance entanglement coefficient y'^ (dashed curve in the insertion); distance 
entanglement coefficient Z (EBI). 



Figure 8: The covariance entanglement coefficient 3^1,2 (87) in the ID "semiresonance" (p = 1) cavity versus 
the "slow time" t, for the Fock (dashed curves) and squeezed vacuum (full curves) initial states of the first 
mode with mean photon numbers vi ~ 1, 50, 1000. 

Figure 9: The mean number of photons in the first and second modes of the ID "semiresonance" (p = 1) 
cavity versus the "slow time" t, for the initial Fock state |1). 



Figur e 10: The purity entanglement coefficient £1,2 (108) (full curves) and distance entanglement coefficient 



Z12 (109) (dashed curves) versus the "slow time" r, for the ID "semiresonance" [p = 1) cavity and the 



initial squeezed vacuum state of the first mode with different mean numbers of photons vi = 1, 50, 1000. 
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